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Using a second-order perturbative Green’s functions approach we determined the normal state 
spectral function A(k, uj) employing a minimal model for ferropnictides. Used before to study 
magnetic fluctuations and superconducting properties, it includes the two effective bands related 
to Fe-3d orbitals proposed by S.Raghu et al. [Phys. Rev. B 77, 220503 (R) (2008)], and local 
intra- and inter-orbital correlations for the effective orbitals. Here, we focus on the normal state 
electronic properties, in particular the temperature and doping dependence of the total density of 
states, A(io), and of A(k, oj) in different Brillouin zone regions, comparing them with existing angle 
resolved photoemission spectroscopy (ARPES) and theoretical results. We obtain an asymmetric 
effect of electron and hole doping, quantitative agreement with the experimental chemical potential 
shifts, as well as spectral weight redistributions near the Fermi level with temperature consistent 
with the available experiments. I 11 addition, we predict a non-trivial dependence of A(cj) with 
temperature, exhibiting clear renormalization effects by correlations. Interestingly, investigating the 
origin of this predicted behaviour by analyzing the evolution with temperature of the k-dependent 
self-energy obtained in our approach, we could identify a number of Brillouin zone points, not 
probed by ARPES yet, where the largest non-trivial effects of temperature on the renormalization 
are predicted for the parent compounds. 

PACS numbers: 71.10.-w, 71.10.Fd, 71.20.-b, 71.20.Gj 


I. INTRODUCTION 

In 2008, the discovery of superconductivity in 
LaFeAsOi-cFa;, with transition temperature T c = 26 K 
at x > 0.05 doping^ prompted great interest in the ex¬ 
perimental and theoretical study of iron-based supercon¬ 
ductors. LaOFeAs-like compounds, denoted the 1111- 
ferropnictide family, served as starting point and were 
followed rapidly by the discovery of similar electronic 
properties in a series of other iron-based families of com¬ 
pounds. In the meantime many new superconductors 
have been found, all of them including quasi-two dimen- 
sional(2D) layers with a square-lattice array of Fe ions^ 
A vast amount of experimental and theoretical effort has 
been devoted to explore this family of compounds, evi¬ 
denced by more than 15,000 articles published since 2008 
according to the latest reviews, see e.g. Refs—£ and ref¬ 
erences therein. Therefore, in the following we will focus 
on the previous research work most relevant to provide 
a proper context for the discussion of the open problem 
which we have addressed, and the new results on the 
normal state spectral properties we obtained. 

It is generally recognised that the ground state pro¬ 
perties of the ferropnictide parent compounds are mostly 
well described by first principles density functional the¬ 
ory (DFT) calculations, indicating that the low energy 
excitations are mainly due to electrons in the Fe-3d 


orbitals For LaOFeAs, calculations were carried out 
almost simultaneously in a number of works—The 
Fermi surface topology and band structure are rather 
similar for the 11 (e.g. FeSe), 111 ( e.g. LiFeAs, 
NaFeAs ), 122 ( e.g. BaFe 2 As 2 , CaFe 2 As 2 , EuFe 2 As 2 ), 
1111 ( e.g. LaOFeAs, LaFePO, SmFeAsO,CeFeAsO) se¬ 
ries of iron-based superconductors— However, as a 
function of doping the agreement between DFT calcula¬ 
tions and ARPES experiments in materials like LaOFeAs 
diminishes,— and unexpected Fermi surface topology 
changes were found as a function of temperature and 
doping— — 

Theoretical predictions showed that the 1111 com¬ 
pounds consist of two quasi-2D Fermi cylinders at the 
zone center (r) and a massive 3D hole pocket at the 
Z-pointi 10 ’ 24 This behavior was observed by ARPES 
and quantum oscillation measurements— On the other 
hand, in the case of 122 compounds, observations re¬ 
lated to the dimensionality of the electronic structure 
revealed quite different behavior, i.e. a more 3D nature 
of the electronic structure was found in electron doped 
BaFe 2 As 2 — »— whereas a quasi-2D electronic structure 
was derived for the case of K-doped BaFe 2 As 2 — Interest¬ 
ingly, Lui et al— showed a transformation of the nature 
of the electronic structure of CaFe 2 As 2 from quasi-2D to 
more 3D as a function of temperature, going from the 
high-temperature tetragonal to the low-temperature or- 
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thorhombic phases. 

Many aspects of the temperature dependence of the 
normal state spectral properties remain largely unstud¬ 
ied in ferropnictides, which motivated our present work. 
Here, we present our study of the changes with tempe¬ 
rature and doping of the normal state electronic struc¬ 
ture, throughout the Brillouin zone (BZ). We com¬ 
pare our results with ARPES experimental and previ¬ 
ous theoretical results, where feasible, and have been 
able to predict interesting non-trivial temperature de¬ 
pendent effects at Brillouin zone points, yet unexplored 
by ARPES. To do this, we employed a minimal mi¬ 
croscopic model, which includes two effective bands^S 
to describe the low energy band structure, as well as 
intra- and inter-orbital local Coulomb interactions, as 
detailed in next section. Similar effective two-orbital 
models have been studied before, focusing on other as¬ 
pects of the problem, mostly on the superconducting pro¬ 
perties. For example, investigations of pairing mecha¬ 
nisms and gap symmetry were reported/** — studies of 
spin fluctuation s 31 ’ 32 ! 34 and a spin-density wave phased 
Mott transition for strong electron correlationslat¬ 
tice and orbital properties*^ - — etc. In our case, we stud¬ 
ied the minimal microscopic two-orbital model using an¬ 
alytical perturbative techniques to take into account the 
effect of electron correlations, in order to determine the 
relevant normal state Green’s functions and the corre¬ 
sponding temperature-dependent electronic spectral den¬ 
sity. We would like to highlight the fact that, in con¬ 
trast to other analytical techniques, in our approach a 
k-dependent self-energy is obtained, also dependent on 
temperature and doping, enabling us to explore normal 
state physical properties throughout the Brillouin zone, 
and thus address new problems. 

Our perturbative treatment for the electronic corre¬ 
lations in ferropnictides is justified by previous estima¬ 
tions of intermediate values for them, both from numer¬ 
ical calculations and experimental results. For example, 
a combination of local density approximation for DFT 
and dynamical mean field theory (DMFT) calculations 
for REOFeAs (RE = La, Ce, Pr, and Nd)*^* estimated a 
Hubbard local intra-orbital correlation value: U ~ 3.69 
eV. Refs,li^i estimated that with U = 4 eV, width of 
Fe-bands Wp e = 3 eV from LDA, and a Hund's coupling 
energy J = 0.7 eV, taking for the inter-orbital correla¬ 
tion V = [7 — 2,7 as in most refs.* 46 ’ 48- — in this case 
yielding V = 2.6eV, the available normal state ARPES 
results could be described, while the same parameters 
were also adopted for studying the superconducting state 
in Ref—. Furthermore, it was concluded^ that U = 4.5 
eV would transform the system into a Mott-insulator, 
in contrast to other LDA+DMFT predictions*^ - — re¬ 
porting that the system remains metallic and does not 
transform into a Mott-insulator even increasing U up to 
5 eV, with Wp e = 4eV and J = 0.7 eV; while for J = 0 
an electronic structure characteristic of much weaker cor¬ 
relations is obtained, with a quasiparticle weight renor¬ 
malization factor Z ~ 0.8— From ARPES, mass renor¬ 


malization factors between 1.5 and 2.5 were estimated in 
122 compounds of the form Sri_ a ,K a; Fe 2 As 2 , depending 
on the Fermi surface sheet *^“—- with similar reports of 
1.3 - 2.1 for 1111 and other 122 compounds ] 25 ’ 58 ’ 59 X- 
ray absorption spectroscopy estimations of U, placed it 
below ~ 4 eV j 60 ’ 61 To describe the smaller experimental 
Fermi-surface areas reported in de Haas- van Alphen ex¬ 
periments in SrFe 2 As 2 r^ a renormalization of the LDA 
band structure was suggested in 2009estimating an in¬ 
terband scattering of magnitude ~ 0.46eV between hole 
and electron bands to explain the reported experimental 
band shifts— Regarding the minimal two-orbital model 
by Raghu et al. of Ref. —, the full bandwidth of the 
non-interacting band structure is about 12 eV, and lo¬ 
cal intra-orbital correlation values U/W = 0.2 — 0.5 have 
been used previously—^ Adding Hund’s coupling J and 
local Coulomb correlations, a Mott transition in the two- 
orbital models was predicted at a critical interaction 
UJW ~ 2.66/(1 + J/U)P- if J/U < 0.01, with a de¬ 
crease of U c for larger J: deducing that the effect of in¬ 
creasing Hund’s coupling on critical U c is similar to the 
effect of having more degenerate orbitals in a multior¬ 
bital Hubbard model*^ effectively weakening the effect 
of correlations^ and stabilizing a metallic state— 

In addition to effective two-orbital models for fer¬ 
ropnictides, other nrultiorbital effective models were 
proposed^ A three-orbital Hamiltonian was constructed 
involving the 3d orbital xz,yz and xy for Fe^, and com¬ 
pared with the two-orbital model of Ref. 39 . Improve¬ 
ment of two shortcomings of the latter were reported, 
related to the relative weights of each orbital on the 
Fermi surface. Also effective four-orbital 41 and five- 
orbital models^ ’ 11 ’ 32 ’ 48 ’ 53 ’ 64 were studied. Nevertheless, 
the effective two-orbital model of Ref. 39 is still recognized 
as a useful minimal model to describe the main features 
of the low-energy physics of ferropnictides— -37 ’ 41 ’ 42 ’ 66 

Our paper is organized as follows. In Section III Al 
we present the microscopic correlated two-orbital model 
adopted for our study of the normal state of Fe-based su¬ 
perconductors, and describe the analytical Green’s func¬ 
tion approach we used to calculate the electronic spectral 
density function and the total density of states including 
the correlations in second-order of perturbations ( fur¬ 
ther details of our analytical calculations appear in Sup¬ 
plementary Appendix A). In Section IlIIl we present and 
discuss electronic structure results obtained at different 
temperature and doping values along the Brillouin zone, 
and compare them with available theoretical and experi¬ 
mental ARPES results for the normal state spectral pro¬ 
perties of ferropnictides. To understand the origin of the 
non-trivial temperature dependence we predict for the 
density of states, in Supplementary Appendix B we an¬ 
alyze the k-dependent electron self-energies obtained in 
our approach, having identified a number of specific Bril¬ 
louin zone points, not probed by ARPES experiments 
yet, where we find the largest non-trivial effects of tem¬ 
perature on the renormalization. We conclude in Sec¬ 
tion El summarizing our main findings with the mini- 
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mal model for pnictides and analytical treatment used, 
prompting for new ARPES experiments which could test 
our predictions, and mentioning possible extensions and 
applications of our work. 


II. MICROSCOPIC MODEL AND ANALYTICAL 
APPROACH. 


A. Correlated effective two-orbital model. 

To describe analitically the properties of ferropnictides, 
we will consider a simplified model which contains the 
minimum number of degrees of freedom preserving the 
essential physics of the problem in 1111 and 122 com¬ 
pounds, as mentioned in previous section. In particu¬ 
lar, a minimal two-orbital model suitable to describe the 
Fermi surface topology of ferropnictides was proposed by 
Raghu et al. in Ref. 30 . The model consists of a two- 
dimensional square lattice, with each site having two 
degenerate orbitals. Raghu et al. in Ref— 2. fitted the 
bare-band tight-binding parameters, to obtain an effec¬ 
tive band structure which, after folding the Fermi surface 
to the two Fe/cell Brillouin zone, exhibits two hole pock¬ 
ets around the F point and two electron pockets around 
the M point, which are characteristic of the majority of 
FeAs compounds—On the large one Fe/cell Brillouin 
zone, which we will be using in this work, the unfolded 
Fermi surface includes a hole-related band, relevant to 
the description of the physics near the T point, and an 
electron-related band to describe the physics near the M 
point. The two-orbital model was also shown to be suit¬ 
able to describe the extended s-wave pairing believed to 
be relevant for ferropnictides, and other details of the 
superconducting phase—^ ~ 37 i 66 

In this work, we will study the doping and temperature 
effects on the electronic structure and spectral proper¬ 
ties, using the following microscopic model for the normal 
state of pnictides: 


U = 'Hv + V int . (1) 

Thus, we will be modelling these compounds by an 
effective extended Hubbard model: consisting of two cor¬ 
related electronic orbitals, described by Hamiltonian Hq 
including the two effective bands proposed by Raghu et 
al. 30 in the large BZ with one Fe atom/unit cell, and local 
intra- and interorbital electronic correlations included in 
Vi n t■ Concretely: 


= J2 [E c (k)clc ka + E d (k)dld k 

k,a 



where the operators c' k and d\. create an electron in the 
respective c and d effective bands, with spin cr 4- and 
crystal momentum k while, following Ref. 3 —, we consider 
for the two bare effective bands: 
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E c ,d$) = e +(k) ± \Je 2 _(k) + el y (k) - /z (3) 

/z denotes the chemical potential: notice that we refer the 
energies to /z, since we will adopt the grand canonical en¬ 
semble in our Green’s function treatment of the system, 
like in Ref. 30 , o denotes the spin degrees of freedom, and 


e± (k) = (4) 

= —4t 4 sin(k x ) sm(k y ) 

e x (k) = —2t\ cos (k x ) — 2 t 2 cos (k y ) — 4f 3 cos (k x ) cos (k y ) 
e y (k ) = —2 1 2 cos (k x ) — 2ti cos (k y ) — 41 3 cos (k x ) cos (k y ) 

where the parameters ti,i = 1,4, denote the hopping 
amplitudes between sites on the square lattice formed by 
the Fe atoms22 in ferropnictides. 

As for the electronic correlations, we assume that 
short-range Hubbard-like electron-electron interactions 
are present, and consider two local correlations between 
the effective orbitals at each site ( as in Refsj 30 ’ 31 ): intra¬ 
orbital Coulomb repulsion U, and inter-orbital repulsion 
V. Therefore V. nt has the following form: 


Vint — 'y ' [U (ni^riii + Ni^Ni^) + V (n ^ + trj.) (N^ + iVq,)] 

i 

( 5 ) 

Here: rii a = c\ a Ci a and Ni a = d\ a d lcr , and i denotes the 
lattice sites. Fourier transforming, I ^ nt can be written in 
fc-space as follows: 

^ int = ~N [ c fcit Cfc 2 t c fc 34 . Cfc i- fe 2 +fc 3 4 . + 

fe1.fc2.fc3 

[4i t c fe2t 4 3t 4 1 -fc 2+ fc3t + 4 lt ^t43t4 1 - 

fc1.fc2.fc3 

+ [4i4. c fc'2l4 3 t^ fc i- fc 2+fc3T + 4i4, c fc24-4 3 4.dfci- 

fcl.fc2.fc3 

( 6 ) 


B. Our electronic structure calculation for 
ferropnictides. 

The total electron spectral density to be compared 
with ARPES experiments is A(k,co), which in our effec¬ 
tive extended Hubbard model is defined by: 

A(k,uj) = A c (k,ui) + A d (k,uj) (7) 

with contributions from both electron bands given by: 


fc2 + fc3t 


fc2 + fc3t 


fc2 + fc3t 
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and in Eq. 1131 


A ccr (fc,w) = ——ImG r a et (k, w); A da (k,uj) = ~-ImF^ et (k 

7 r 7r 

( 8 ) 

Here: G r a et {k,uj) and F^ et {k, w) denote the retarded 
Green’s functions corresponding the c and d electrons in 
our two effective bands, i.e. 

G r a et (fc, w) = G(k,u) + iS) =< Cko ; c\ a > (w + id) (9) 


,) r 4 (fc 1 ,/c 2 , k,u) = r££ w (k 1 ,k 2 ,k,uj) 

r 5 (k 1 ,k 2 ,k,u) = r^ iCT (fci, k 2 ,k,u) 

r 6 (fci,fc 2 ,fc,w) = r^ 5 :(fci, k 2 ,k, to) (15) 

having introduced the following notation, above, for sim¬ 
plicity: 


F£ et (k,u}) = F(k,oj + i5) =<d fco .;4 <T > (w + *5) (10) 

where <5 is an infinitesimal positive number. Integrating 
the spectral density over the 1st Brillouin zone, yields the 
total density of states (TDOS), or local spectral function: 

A(uj)= Y A(k,u) ( 11 ) 

kcBZ 

We calculated the temperature-dependent Green’s 
functions using Zubarev’s equations of motion (EOM) 
formalism ) 67 : 68 i.e.: 

u) < A; B »= —<{H,H})+ < [A, H\;B^> (w) 

Z7T 

where A and B are fermionic operators, < > is 

the time Fourier-transform of the retarded Green’s func¬ 
tion —id(t — t')(A(t)B(t') + B(t')A(t)), where the time- 
dependent operators appear in Heisenberg representa¬ 
tion, and the expectation value is calculated using the 
appropriate statistical ensemble at finite temperature T 
(or the ground state of the system, at T=0 )££. In our 
case, we study the normal state of the system at tempe¬ 
rature T, and in particular evaluate the required expec¬ 
tation values in the paramagnetic normal phase, using 
the grand canonical ensemble. 

We obtained the following exact coupled set of equa¬ 
tions of motion for G a (k,u}) and F a (k,co), respectively: 


w 


-E c (k)]G ir (k,w) = — + Y 


ki,k2 L 


^V 1 (k 1 ,k 2 ,k,uj) 


V V 

+—T 2 (k 1 ,k 2 ,k,u}) + —T 3 (ki,k 2 ,k,u) 


( 12 ) 


[w - E d (k)\ F a (k, w) = ^ + Y 


ki,k2 


Ti(ki,k 2 ,k,u) 


V V 

+—r 5 (k 1 ,k 2 ,k,u) + —r 6 (fci,fc 2 ,fc,w) 


(13) 


where in Eq. [T2]we denoted: 


ri(/ci,fc 2 ,fc,w) = r^ w (ki,k 2 ,k,u) = 

= Ck 2 ,ac\. l W Ck 1 —k2+k,W: c k cr ^ ( w ) 

T 2 (k 1 ,k 2 ,k,co)=T^ a (k 1 ,k 2 ,k,u J ) 

r 3 (fcr ,k 2 , k, w) = T l%(k u k 2 , k, u) ( 14 ) 


T''0!,/3,'7 
& a. &/3 


(fcl,fc 2 ,fc,w) EE < Co 


^fcl„ 


c lk 1 -k 2 + k,v 1 i C k<j '* E> 


(HI) 

where c a ,c^,c 7 describe c- or d-annihilation operators 
and cx a ,<Tg,cr 7 their spins (a, or a), as required by each 
equation (more details in supplementary Appendix A). 

Notice that the Hartree-Fock solution to this problem 
would be obtained by closing this set of equations of mo¬ 
tion in first-order, with a mean-field approximation of Tj 
( i = 1,6) in terms of G a (k,u>) and F a (k,u}), as detailed 
in supplementary Appendix A. 

In our work, instead, we calculated the equations of 
motion for the three I\ (* = 1,3) functions coupled to 
G(fc,w) in first order as stated in the previous equa¬ 
tions; and proceeded likewise for Tj (i = 4,6) coupled 
to F(fc, w). Each of these second-order equations of mo¬ 
tion introduces three new coupled higher-order Green’s 
functions in the problem. 

To close and solve the coupled set of equations of mo¬ 
tion in second-order, we used the following approxima¬ 
tion: all higher-order Green’s functions introduced in 
each subset of equations of motion for Tj (i = 1,3), 
were approximated in mean-field in terms of T,; (i = 1, 3) 


and G(k, to) (introducing proper average values), yielding 
a 4 x 4 closed set of equations of motion ( for details, 
see supplementary Appendix A). We proceeded likewise 
for the subset of equations for Ti (* = 4,6) related to 
F{k,uj). Notice that this RPA-likc approximation leads 
to a second-order closed system of eqs. of motion for G, 
lb, T 2 , 1?3; and an analogous one for F, lb, Ts, 1^; which 
after lengthy calculations we could solve. In supplemen¬ 
tary Appendix A we detail the full expressions obtained 
for these Green’s functions, which might be used to calcu¬ 
late other normal state properties, and in supplementary 
Appendix B we discuss the corresponding k-dependent 
electron self-energies, which are obtained in our level of 
approximation. 

The decoupling scheme proposed for the higher order 
Green’s functions, which appear coupled in the EOM, 
accounts for the self-energy effects of the single-particle 
Green’s function. These self-energy effects may describe 
both intra- and inter-orbital coupling with equal and 
opposite spins at lowest order (charge and spin fluctu¬ 
ations), thus leading to band splitting effects, which are 
not a mere renormalization of the energy spectrum. This 
fact is well known for other Hubbard-like models (see e.g. 
Refs j 68 ’ 69 ), and the band splitting effects may lead to a 
metal-insulator transition or to charge and spin instabi¬ 
lities. 
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We shall be interested in describing the behavior of 
this many-particle system at finite temperatures. For 
a system in thermodynamic equilibrium the expectation 
value of any operator may be computed by using the 
grand-canonical ensemble. To complete our solution in 
RPA approximation, the total electron band filling n (or 
the chemical potential n) at temperature T is determined 
self-consistently by the following equation: 


n(n) 


dk dco 

(27r) 2 27T 


A c (k,i 


1 




+A(z(/c, U>) 


1 


(17) 


Thus, a major self-consistent set of coupled equations 
has to be solved, since the total spectral density function, 
which we obtain by solving two coupled sets of Green’s 
functions, is required and will be integrated over the Bril- 
louin zone (fc) and all energies w, weighed by the Fermi 
occupation factors of each effective band at temperature 
T . The thermodynamic state of the system is defined 
by the parameter /i and jd, the inverse temperature mea¬ 
sured in energy units, i.e. fd = where ks is Boltz¬ 
mann’s constant. Zero temperature, or /3 —>■ oo, describes 
the ground state of the system. Notice that in this model, 
half-filling (n = 2) corresponds to the ferropnictide par¬ 
ent compounds. 

As detailed in supplementary Appendix A, to evaluate 
the Green’s functions G CT (fc,w) and F a (k,oj), we need to 
perform double and triple summations over the first Bril- 
louin zone of the crystal lattice. For simplicity, we have 
assumed a square lattice, as done previously^ and have 
performed the BZ summations using the Chadi-Colien 
BZ sampling method. 70 


III. RESULTS AND DISCUSSION 

Below, we present electronic structure results at diffe¬ 
rent temperature and doping values, which we obtained 
for the normal state of ferropnictides with the extended 
Hubbard model and our analytical treatment presented 
in previous section, and compare them with available the¬ 
oretical and experimental ARPES results in ferropnic¬ 
tides. The self-consistent sets of equations obtained for 
the Green’s functions detailed in previous section, were 
solved numerically using the following tight-binding in¬ 
teraction parameters: t\ = —1 eV, t .2 = 1-3 eV , = 

t 4 = —0.85 eV, adopted in Refill! for the two effective 
bands in the large one Fe/cell Brillouin zone: BZ which 
was also used throughout our work. In the figure captions 
we specify the precision used in the Cliadi-Cohen sam¬ 
pling method 70 for Brillouin zone summations, including 
the order v used. We found good convergence in most re¬ 
sults using the seventh order {v = 7) of the Chadi-Cohen 
method, 70 but the spectral function and total density of 
states results here presented were obtained using ninth or 
eighth order, respectively, for improved accuracy. Note 
that v = 8-th order for the square-lattice BZ sampling 


implies using 8,256 special BZ points, while v = 9 implies 
using 32,896 ones^S 

We will start by analizing in subsection IIII Al the effect 
of correlations in our model, comparing our calculations 
of the density of states to previous known results, like 
experimental data now available for 1111 and 122 com¬ 
pounds, in order to obtain a qualitative estimation of the 
relevant magnitudes of U and V in the model. 

Next, in subsection IIII Bl we exhibit the momentum 
dependence of the spectral density function along sym- 
mfi'.ry paths of the square lattice Brillouin zone (BZ), 
Hcdmparing our results with reported ARPES data. No¬ 
tice that in order to allow for direct comparison with the 
published ARPES intensity results, the same treatment 
is adopted here for the data: i.e. all spectral density data 


here presented ( denoted A(k,co), here) have been con¬ 
voluted with an energy resolution of 30 meV, and also 
multiplied by the Fermi-Dirac function /fd(w) at tem¬ 
perature T as in the experimental references cited. 

Subsequently, in subsection IIII Cl we study the effects 
of hole and of electron doping in terms of our model, and 
compare them with the electronic properties of ferrop¬ 
nictides observed in ARPES. Finally, in subsection IIII Dl 
we study and compare with available experimental data 
the temperature dependence of the total density of states 
and the spectral density at the most relevant BZ points. 


A. Effect of intra- and inter-orbital Coulomb 
interactions on the total density of states A(uj) 

As discussed in the Introduction, previous theoreti¬ 
cal and experimental research works estimate interme¬ 
diate values for the electronic correlations, i.e. most 
of them placing the local intra-orbital Hubbard repul¬ 
sion in the range U ~ 3.69 eV for ReFeAsO (Re= Ce, 
Pr, Nd)^ to U = 4.5 eV for LaOi_ 3 ;F x FeAs , 47 i 52 while 
X-ray absorption experiments suggest U < 4 eV for 
LaOi-^F^FeAs^i and BaFe 2 As 2 ^£ For the inter-orbital 
correlation V similar values are estimated in calculations 
for Fe-Se compounds^ 

We analized the separate effects of the intra- and inter¬ 
orbital correlation parameters included in the model. 
First, the total density of states A(u), is shown in Fig. 
Oja), for different values of intra-orbital U and a fixed 
inter-orbital: V = 3.5eU, at half-filling n = 2 and tem¬ 
perature T = 20K. We find a large spectral weight re¬ 
duction at the Fermi level in the range U = 3.7 — 3.9 eV, 
with an important spectral weight redistribution around 
the Fermi level. 

In Fig. [D)b), we focus on the effect of the inter-orbital 
Coulomb interaction V on the density of states. Note 
that the TDOS is almost independent of V, in the range 
of U values analized. We thus confirm that the most 
relevant Coulomb correlation in 1111 and 122 systems is 
in fact the local intra-orbital repulsion U, in agreement 
with previous research work. 7 ' 3 ' 

In the following, we will focus on the dependence of the 















B Momentum dependence of the spectral density A(k , w) 
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FIG. 1. a) Dependence of total DOS A(lo) on the intra-orbital 
repulsion U, for fixed V = 3.5eV ; b) Effect of the inter-orbital 
local Coulomb interaction V on the total DOS, for fixed U = 
4.5eF. At half-filling: n = 2, and temperature: T = 20 K; 
and as in Ref.—: ti = —1.0, t 2 = 1.3, f ,3 = I 4 = — 0.85 (in 
eV). Chadi-Cohen BZ summations order: v = 8 . Note that 
oj is measured w.r. to p, in our treatment. All energies are 
given in eV, henceforth. 


electronic structure on crystal momentum, doping and 
temperature, fixing the local Coulomb interactions as: 
U = V ~ 3.50 eV, typical intermediate correlation val¬ 
ues characteristic of Fe-pnictide compounds as previously 
discussedi 52 ! 60 ’ 61 : 71 Being the total bandwidth W ~ 12 
eV for the non-interacting effective band structure of 
Ref—, the value U = 3.5 eV corresponds approximately 
to 0.29 W. 


B. Momentum dependence of the spectral density 

A(k, lo) 

Our analysis of the two-orbital extended Hubbard 
model included the calculation of the one-particle spec¬ 
tral function. In this section, we will show our electronic 
structure results along two BZ paths, namely along T - 
X in Fig. (2](a), and along T - M in Fig. [IJb), where: 
T = (0,0), X = ( 7 r, 0) and M = (7r, 7r). These paths were 


chosen as they allow us to compare our results with avail¬ 
able ARPES data, in particular those of Refs j 65 i 72 ~ — . 
Concretely, we here show the evolution of the spectral 
function A{k , to) as a function of momentum k , at low 
temperature: T = 20 K, and at half-filling: n = 2 (i.e. 
for the parent compounds). 



-0,6 -0,4 0,2 0 


(o (eV) 



FIG. 2. a) Spectral density A(k, lo) (i.e. convoluted with 
an exp. resolution of 30 meV, and multiplied by /fd(w), as 
in expts.) as a function of energy lo, for different BZ points 
(shown vertically displaced) along the symmetry path from T 
(bottom curve) to X (top curve), b) A(k,uj) as a function of 
U, for different BZ points along T-M. At T = 20 K, and n = 2 
(as in Ref.—). Here: U = V = 3.5eF; p = 0.15eV. Chadi- 
Cohen BZ summations order: v = 9. Other parameters as in 

Fig.m 

In particular, in Fig. Ela) notice that a relatively flat 
multiple-peak structure near T, evolves into a sharp peak 
near the Fermi level at X. In accordance with ARPES 
experiments along r — X, e.g. in LaOFe(P,As),(Fig. 
4 of Ref—) in K a; Fe 2 -j / Se 2 ( Fig. 2 of Ref—) and 
BaFe 1 . 7 Coo. 3 As 2 compounds,( Fig. 2 of Ref—) during 
this evolution, the main peak describes an S-like trajec¬ 
tory including an intermediate regime with a two peak¬ 
like structure, though in experiments the width and in¬ 
tensity of these features are material-dependent. 

Along r - M, in Fig. [2])b) notice how the broad rel¬ 
atively flat structure near T and near M, for intermedi¬ 
ate BZ points develops a peak-like structure which de¬ 
scribes a trajectory similar to that reported in ARPES 
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experiments^ concretely, Fig. 1 of Ref.— and Fig. 1 of 
RefJ£) in particular, it resembles the evolution of the 
ARPES peak in LaFeAsO assigned to two bulk-related 
energy bands in RefJ^, where a second peak-like struc¬ 
ture closer to the Fermi level but related to surface-bands 
is also shown. Very recently, ARPES in FeSe single crys¬ 
tals reported^ very similar energy distribution curves 
along F - M in their Figures 1 and 2. 

Notice that the spectral density results shown in Figs. 
Ha) and[2jb) include renormalization effects due to the 
correlations which are taken into account by our analyti¬ 
cal approach, concretely by the electron self-energies ob¬ 
tained, discussed in Supplementary Appendix B, which 
depend explicitly on crystal momentum k , as well as tem¬ 
perature and doping. Band renormalization effects such 
as these, implying changes to the non-interacting electron 
band structure of the minimal model proposed in Ref. 30 
(in particular, to their Fig. 2(a)), have also been reported 
in several ARPES experiments: e.g. in Refs . 20 ’ 54 ! 74 . 


C. Effect of doping on the electronic structure 

We will now refer to the effects of doping on the total 
density of states, and in particular the degree of accuracy 
presented by the results obtained with our approach, tak¬ 
ing advantage of the experimental data available for com¬ 
parison, such as Refs. 76- ' 7 . In Supplementary Appendix 
C we complement these results, showing also the effects 
of doping on the spectral density function at the T and 
M symmetry points. We carried out a systematic study 
covering a broad range of hole and electron doping values, 
which respectively correspond to less or more electrons 
in the system than present at half-filling (n = 2), which 
represents the parent compound. In our calculations we 
fixed the band filling n, while the chemical potential p 
was obtained self-consistently according to Eq. 17. 

In Fig[3] (a) we exhibit a comparison of the effect of 
hole vs. electron doping on the renormalized total den¬ 
sity of states. In agreement with ARPES experiments in 
Ref.' 6 (in particular, their Fig. 2(e)), we find that upon 
increasing hole-doping a rigid-band-like shift of all peaks 
towards the Fermi level ensues, thus increasing spectral 
weight there ( for more details, see Suppl.Appendix C. 
In particular, in Figs. C.5 and C.6 the spectral density 
results at M are shown, and their comparison with ex¬ 
periments is discussed in detail). Meanwhile, when the 
system is doped with electrons, n = 2.12 and n = 2.25 
in the figure, the bottom of the electronic structure is 
seen to be shifted to lower binding energies. Notice that 
for higher electron-doping values, e.g. n = 2.25, a sharp 
intense peak appears at the Fermi level. 

For comparison with the results of Figure 3 of Ref. —, 
in our Fig [3] (b) we plot the chemical potential shift: 
Sp = p(n) — p(n = 2) with respect to the undoped 
compound obtained in our approach as a function of 
the carrier (dopant) density: 5n = n — 2. As indi¬ 
cated in Ref. —, experimental chemical potential shifts, 



FIG. 3. Electron vs. hole doping, a) Effect on the total DOS. 
b) Chemical potential shift as a function of carrier (dopant) 
density, with respect to the parent compound. For the exp. 
values shown, obtained from reported^ core-level shift data 
ApcLS (solid triangles) and valence-band shifts A pvbs (solid 
circles), the carrier density is given as x/2 (per Fe): either for 
BaFe 2 -r E Co a: As 2 or Bai_ a; K a; Fe 2 As 2 , respectively for electron 
or hole doping. At T = 20 K, U — V = 3.5eU. Other 
parameters as in Fig. |T] 


which we denote A pcls, can be obtained from the do¬ 
ping dependence of the measured As-3d core level shifts 
( A Ecls ) reported in their Fig.3.b: by using the for¬ 
mula A pcls = -AEcls + (A Vm + AEr) , where 
( A Vm + AEr) accounts for doping-related changes in 
the Madelung potential and in the core-carrier screening, 
respectively, assuming that changes of As-valency upon 
doping are negligible. Using the As-3d core level shifts 
and the differences (A Vm + A Er) reported in Fig.3(b) 
of Refj—, we calculated the experimental chemical po¬ 
tential shifts derived from core-level shift data A pcls 
corresponding to the seven compounds experimentally 
studied, and have plotted them in our Figure HR)- A 
comparison with the theoretical results of our approach, 
shows that the asymmetry and especially the different 
signs exhibited by the chemical potential shifts for hole 
and electron doping are described, in agreement with 
experiments—2’ 76 Notice that our theoretical results also 
describe the correct monotonous dependence of the chem¬ 
ical potential shift as a function of doping for all samples, 
and quantitative agreement is also obtained within the 
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experimental error bars of A pcls for all the electron- 
doped samples as well as for the sample with the high¬ 
est hole-doping. For completeness, in Figure EJb) we 
also include the experimental chemical shifts, denoted 
by A pvbs , obtained as A pvbs = —A Evbs from the 
valence-band shift data (A Evbs) reported in Fig.3(c) 
of Ref—, which are consistent with A pcls obtained 
from the measured As-3 d core level shifts. The reported 
electron-hole asymmetry in Ba-122 pnictides/ 6 in parti¬ 
cular in BaFe 2 As 2 , has been also observed in the super¬ 
conducting phase. 20 ' 76 ’ 7 9'8£ We checked that the chemical 
potential shift is almost independent of V, consistently 
with our DOS results discussed in Section UlIAl 

Let us first stress one result: notice that in spite of the 
band renormalization effects due to correlations known 
to be included in our level of approximation ( as shown 
in our self-energy analysis of Suppl.Appendix B), we ob¬ 
tain electron-hole asymmetric effects such as experimen¬ 
tally observe d 76 ' 77 i.e. the bottom of the band is shifted 
to lower, or higher, binding energies with respect to the 
Fermi level upon, respectively, electron or hole doping, 
as one would expect if the electronic structure is not 
changed by doping. 

Next, we would like to stress that, even though one 
might expect the values of the chemical potential shifts 
upon doping to be material-dependent, we find that our 
results in Figj3](b) agree better with the reported data 
on Ba-122 compounds (Figure 3 of Ref. 7 —) than the LDA 
calculations mentioned in Ref.' 6 : which required divid¬ 
ing their calculated chemical potential shift values by a 
factor of 4, in order to achieve quantitative agreement 
with experiments. Previously, theoretical estimations^ 
of an LDA band structure renormalization, amounting 
to a factor of 2 division, had been indicated as necessary 
to describe with LDA the observed energy shifts of the 
electron and hole bands— 

Thus, our analytical treatment for the effective two- 
band model yields electronic structure results describing 
not only the observed opposite chemical potential shifts 
for electron and hole-doping, but the chemical potential 
shift values we obtain are in quantitative agreement with 
the experimental results of Fig. 3 of Ref. 7 -, respectively 
for all electron doped samples in BaFe 2 - x Co x As 2 and for 
the sample with highest hole doping in Bai_ x K x Fe 2 As 2 


As mentioned in the Introduction, Lifshitz transitions 
have been reported for ferropnictides, in particular as 
a function of doping— i- 22 ^ We have checked that the 
correlated two-orbital model indeed undergoes changes 
of topology of the Fermi surface upon doping increase, as 
exhibited in Figure C.8 of Suppl. Appendix C. A detailed 
study of this topic is out of the scope of the present work, 
though we plan to address it in the future. 


D. Effects of temperature on the spectral 
properties 

Our analytical approach is specifically well adapted to 
shed light on the temperature dependence of the normal 
state spectral properties throughout the Brillouin zone, 
subject largely unstudied in ferropnictides, on which we 
were able to find interesting predictions. 



FIG. 4. Temperature dependence of A(u). Parameters: U 
= V = 3.5 eV , n = 2. Other parameters as in Fig. [I] Inset: 
partial c- and d-band densities of states at 40 K. 

First, we will focus on the effect of the temperature 
on the total density of states. We found only low tem¬ 
perature density of states ARPES data reported for the 
parent compounds, e.g in Refs— >22, to compare the re¬ 
sults of this section, while at room temperature we only 
found a hard X-ray photoemission DOS result reported 
for undoped Ba-122.— To allow a direct comparison with 
experimental data, in this section all density of states 
results exhibited, denoted A(w), have been convoluted 
with an energy resolution of 30 meV and multiplied by 
/ fd ( w), as in experiments.— ! 81 ! 82 . 

In Figure[4]we show the total DOS we obtain for several 
temperatures varying from 10 K to 100 K. As an inset, 
we plotted the corresponding partial densities of states at 
T = 40 K, corresponding to the two effective TB bands22 
used in our model, denoted by c and d and defined in Eq. 

[HI 

In Fig. [4] notice that with our simplified corre¬ 
lated model for pnictides the total DOS near the Fermi 
level consists of a sharp intense peak, which the in¬ 
set and our self-energy results presented in Supplemen¬ 
tary Appendix B.l show is mostly due to the renor¬ 
malized c-effective band, i.e. obtained in our approach 
by renormalization of the non-interacting c-band from 
Eq. [3] proposed by Raghu et al—, which is the low¬ 
est energy effective band extending from -8 to 4 eV. 
Our result agrees with previous work, which indicates 
that the peak originates mainly from Fe-3d states, in 
particular: ARPES experiments at low temperature in 
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LaOFeP and LaOFeAs compounds (Figure 1 of Ref—) 
at T = 20 K interpreted with DFT+LDA calculations, 
ARPES at T = 28.7K in FeSei_ x (Figure 2 of Ref^), 
and hard X-ray photoemission experiments in BaFe 2 As 2 
(Figure 2 of Ref—) at 20K and 300K. The low tem¬ 
perature results agree also with LDA+DFMT calcula¬ 
tions in SmOi_ x F x FeAs,— LDA+DFT calculations in 
LaFeAsOi_a,F x? iS BaFe 2 As 2 ^ and LiFeAs— 

Regarding the evolution with temperature shown in 
Fig. 01 our calculations predict that this intense sharp 
peak near the Fermi level should be clearly observed for 
temperatures T < 40K. By increasing temperature, we 
find that the peak near the Fermi level is shifted towards 
lower binding energy, diminishes intensity and becomes 
broader, eventually developing additional structure. This 
agrees with the observations in BaFe 2 As 2 (Figure 2 ofSi) 
where, from the sharp peak near the Fermi level at 20 
K, at room-temperature a broad hump located at lower 
binding energies was reported. 

Notice in Fig. 0 ] that we also find new temperature 
dependent peaks in the DOS due to correlation effects, 
as our self-energy analysis in Supplementary Appendix 
B confirms. Concretely, if one compares the density of 
states curves corresponding to various temperatures in 
Fig. 01 at 10 K the DOS is dominated by a single peak 
near the Fermi level, at 40 K we obtain two prominent 
peaks, shifted further below the Fermi level w.r. the T= 
10K result, and three sizeable peaks of the density of 
states ( shifted still further below the Fermi level) are 
present at T= 100 K. Therefore, we are predicting a non¬ 
trivial temperature dependence of the total density of 
states due to correlations, which it would be interesting 
to investigate experimentally. 

Furthermore, in connection with the location of the 
peak near the Fermi level, which may be material- 
dependent: in Fig. 0]we find it centered at ui ~ —10 meV 
at 30 K, which is similar to the peak position estimated 
from fits in Ref.— in FeSei_ x , placing it at ui = —15 meV 
at 30 K, in particular their Figure 2(d). 

Next, we will analize the evolution with temperature 
of the spectral density, analyzing the changes in A(k,ui). 
In first term, we will fix k at the relevant high-symmetry 
points of the BZ which have been studied in ARPES 
experiments j 15 i 29 

In Figure 0](a), we show the evolution of the electronic 
structure with temperature at T. Notice that a num¬ 
ber of prominent peaks emerge as temperature decreases. 
This is qualitatively what was observed in ARPES exper¬ 
iments in LaFeAsO (Figure 7(a) of Ref—) and CeFeAsO 
(Figure 4 of Ref. —) compounds, where three peaks were 
also reported in the same energy range, though our re¬ 
sults differ in the relative intensities. Our results also 
agree qualitatively with recent ARPES experiments at 
T closer to the Fermi level in NaFe 0 . 95 Co 0 . 05 As (Figure 
3(a) of Ref.—) and Cai_ rz .Na x Fe 2 As 2 (Figure 2 of Ref.—) 
compounds, who focus on the spectral density reduction 
at the Fermi level upon temperature increase. 

In Figure 0][b) we show the temperature evolution of 




FIG. 5. a) Temperature dependence of A(k, w), at F; b) Tem¬ 
perature dependence of A(k, u>), at M. Same energy range 
used for the ARPES data of Ref.—. Parameters: U = V = 
3.5 , n = 2. v = 9. Other parameters as in Fig. [T] 


the spectral density at M, to be compared with the re¬ 
ported ARPES data in Figure 3 of Ref— and Figure 
4(b) of Ref—, who focus on the energy range [—0.2, 0]eV. 
The description of the low temperature ARPES data at 
M with the effective two-orbital model is not as good as 
for r. Though the spectra at M in the energy range 
[—0.1,0]eV exhibit a two-peak structure at low tempe¬ 
rature which evolves to a broad single peak at higher 
temperature, in agreement with ARPES experiments in 
CeFeAsO (Figure 3 of Ref—) and EuFe 2 As 2 (Figure 4(b) 
of Ref. 88 ) compounds, we find differences (e.g. a number 
of additional peaks at low temperature below -0.1 eV) if 
we compare our results with experiments in a wider en¬ 
ergy range. The temperature evolution of the electronic 
structure at T exhibits less relevant changes than at M, 
in agreement with experiments. 

In any case, the spectral function results at the high 
symmetry BZ points probed by ARPES, which are pre¬ 
sented in Figs. 0](a) and[5][b), mainly exhibit the usual 
effects of thermal broadening, in correspondence with the 
negligible temperature-dependence of the self-energy at 
these points, discussed in Supplementary Appendix B. 

Nevertheless, as discussed in connection with Fig. 0] 
our model predicts a non-trivial temperature-dependence 
of the total density of states, to which the whole BZ con¬ 
tributes. To understand this fact, we explored through- 









D Effects of temperature on the spectral properties 


10 



b) k 2 = (0.6, 0.6)rc 



co ( eV ) 



cific BZ points ( k\ = ( 0 . 557 T, 0 . 557 r), k 2 = ( 0 . 67 T, 0.67 t) 
and ks = ( 0 . 567 T, 0 . 3 l 7 r)), which we could identify as ex¬ 
hibiting the largest non-trivial temperature dependent 
renormalization effects due to correlations (for details, 
see Supplementary Appendix B). Notice also that some 
non-Fermi liquid features are visible at specific BZ points 
at higher temperatures: in fact, Figure[ 6 ](b) shows that at 
T = 100-A the spectral density function at &2 is mostly 
dominated by the incoherent contribution, a strong re¬ 
duction of the quasiparticle peak having taken place. 
Our prediction could be tested by temperature dependent 
ARPES experiments, and these effects should also be 
displayed by correlated multi-orbital models with more 
effective bands (though perhaps at larger correlation val¬ 
ues). In Supplementary Appendix C.l, we discuss how 
doping affects the results discussed above for n = 2 , and 
in particular how it reduces the non-trivial temperature 
dependent renormalization effects by correlations in the 
density of states, with differences between electron and 
hole doping. We also show how doping affects the k- 
dependence discussed for the parent compounds in con¬ 
nection with Fig. EB Our results indicate that a a sub¬ 
tle interplay between correlations, temperature and do¬ 
ping gives rise to the non-trivial temperature dependence 
effects in the spectral properties. 

Finally, the quasiparticle weights Z(k) at the Fermi 
level, defined as Z - 1 (/c) = 1 — -^ReE c (k, w)| w= o, are 
shown at different BZ points in Fig(3 exhibiting the 
non-trivial effect of temperature along the BZ: clearly 
larger at some points like k$ = ( 0 . 567 T, 0 . 3 l 7 r), which 
would be interesting to probe by ARPES. Note that the 
mass renormalization (-^ oc Z) obtained in our appro¬ 
ach corresponds to a correlated metal, while the range 
of quasiparticle weights obtained along the BZ agrees 
with LDA+DMFT results^ who estimated a quasipar¬ 
ticle weight renormalization factor Z ~ 0.8 for J = 0, 
as well as with T=0 slave-spin mean-field studies of two- 
and four-orbital models.— 


FIG. 6. Temperature dependence of A(k,cj), at the BZ points 
respectively indicated in each plot. Parameters: U = V = 3.5 
eV , n c = 1.45 and rid = 0.55 ( n ■** 2), v = 9. Other 
parameters as in Fig. |T| 


out the whole BZ the dependence on temperature of 
the k-dependent self-energy obtained in our approach, 
and in Supplementary Appendix B exhibit our results. 
Interestingly, for the parent compounds (n = 2) we 
have been able to identify a number of specific BZ 
points, not yet probed by ARPES, which exhibit size¬ 
able and non-trivial temperature dependent renormal¬ 
ization effects, allowing to explain the temperature de¬ 
pendence predicted for the total DOS. In particular, in 
Figure [G] we show the evolution of the spectral den¬ 
sity function with temperature predicted at three spe- 


IV. CONCLUSIONS 

Though a large amount of experimental and theoretical 
research work has been done since 2008, the true nature 
of electron correlations in ferropnictides remains an open 
question. Here, using a simplified model based on two 
correlated effective bands and an analytical approxima¬ 
tion to decouple the equations of motion for the elec¬ 
tron Green’s functions, we determined the normal state 
spectral density function and the total density of states. 
With model parameters in the range relevant for iron- 
based superconductors, we could describe: (i) the mo¬ 
mentum dependence measured in ARPES along r — X, 
and the main features along T — M; (ii) an asymmetric 
effect of electron and hole doping, with the character¬ 
istics observed in ARPES experiments. In particular, 
the band structure renormalization due to the correla¬ 
tion effects captured by our analytical approach, allows 
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FIG. 7. Quasiparticle weights at the Fermi level, Z(k), shown 
at a set BZ points (located as indicated in the inset), for 
temperatures: T = 40, 60,100 K. U = V = 3.5 eV and other 
parameters as in Fig(T] ko = (7t/ 2, 7r/2), fci = (0.557T, 0.557r), 
&2 = (O. 671 ", 0.67r), k$ = (0.567T, 0.3l7r),fc4 = (jv, 0.57r) and ks = 
(0.5tt,0). 


to obtain a direct agreement of our calculations with the 
experimental chemical potential shifts reported, with the 
correct order of magnitude in contrast to LDA calcula¬ 
tions. (iii) We found that an important spectral weight 
reduction near the Fermi level takes place if temperature 
is increased, and the evolution with temperature of the 
total density of states which we predict is consistent with 
the relatively few available data from ARPES and hard 
X-ray photoemission experiments, (iv) The conduction 
electron mass renormalization obtained corresponds to a 
correlated metal, while the range of quasiparticle weights 
at the Fermi level which we obtain along the Brillouin 
zone agrees with theoretical LDA+DMFT results which 
estimated Z ~ 0.8 in the absence of Hund coupling. 

Among the results obtained, we want to stress an inter¬ 
esting prediction which might be tested experimentally: 
our work not only predicts a non-trivial temperature de¬ 
pendence of the total density of states, due to correlation 
effects. We could trace the origin of this behaviour to 
specific Brillouin zone regions where the effect of tem¬ 
perature on the renormalization is amplified, by explor¬ 
ing the k-dependent self-energy obtained in our approach 
throughout the Brillouin zone. It would be interesting to 
have this prediction tested, by temperature dependent 
ARPES experiments in the specific BZ points we explic¬ 
itly identified for undoped compounds. 

Our calculated Green’s functions might be also em¬ 
ployed to evaluate other normal state physical properties 
of iron pnictides, at our level of approximation: which 
has the advantage of enabling to describe temperature 
and k-dependent results, in contrast to other theoretical 
techniques where the self-energy is local. 

We have shown that our results describe well most of 
the main common features of the bulk electronic struc¬ 
ture observed, not only 1111 and 122 ferropnictides, but 
also in FeSe compounds which we believe is related to the 


fact that the two correlated effective bands of the model 
capture the relevant details of the 3d-Fe orbitals near the 
Fermi level, lying the As or Se orbitals further below in 
energy. 

The simplified pnictides model and analytical approach 
used in the present work, are suitable for extensions to 
describe details of magnetotransport in FeSe compounds, 
as well as adapting them to describe electronic properties 
of BiS 2 -based layered compounds. Our approach could 
also be useful to address the nature of correlations in 
Hund metals and its relationship with Mott-physics in 
other multi-orbital systems^, problems we plan to study 
in the future. 
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Appendix A: Analytical approximation: decoupling 
of the Green’s function equations of motion. 

As mentioned in section III B1 for the model Hamil¬ 
tonian given by Eqs. ??, [2] and [5] we calculated the 
equations of motion (EOM) of the Green’s functions 
G(k, uj) and F(k, to), respectively corresponding to c- and 
d-electrons defined in Eqs. [9] and [TUI using Zubarev’s 
formalism^!. We obtained the set of exact equations of 
motion detailed in Eqs. [T2] and QjJ] Furthermore, the 
equations of motion for the three Fj (i=l,3) functions 
coupled to G(k, w), and for Tj (i=4,6) coupled to F(k, w), 
introduce coupling to new higher-order Green’s functions 
in the problem. In this appendix we will detail the RPA- 
like analytical approximation used to obtain the Green’s 
functions G(k,u) and F(k,ui), required to calculate the 
spectral density and total density of states defined in Eqs. 
[3 and El 












1 First order solution of EOM: Hartree-Fock decoupling. 


12 


1. First order solution of EOM: Hartree-Fock 
decoupling. 

First, we will explain how one could close and solve the 
coupled set of equations of motion in lowest order, which 
would correspond to the Hartree-Fock approximation. 

The approximation to be used in Eqs. [12] and [13] for 
the Fi (i = 1,6) Green’s functions which appear coupled 
to G(k,oj) and F(k, w), respectively, basically consists in 
approximating F* ( i = 1,6) as proportional to G(fc,w) 
and F{k,ui ), with mean values of the occupation num¬ 
bers of the c- and d-bands appearing as proportionality 
factors, thus obtaining a closed system of equations of 
motion. Below, we show the approximation explicitly for 
T 1 (ki,k 2 ,k,uj): 


Fi(fci, k 2 , k,u>) = r“£-(fci,fc 2 ,fc,w) = 

=< C C/c 2 ,<T C fci c k<7 '^ > ( W ) 

^k\ ,<7 Gci—fc 2 +fc,<f * >< ^. C/c 2 ,(Ti ^ (^0 
Ck 1 &Ck 1 —k2+k,W > &ki ,fci— fc 2 +fc ^ Ck 2 ,cr\Ckcr ^ ( W ) 
n - J< C Tlk i^ > $k,k2 (Al) 


= 4,6) one has: 

r 4 (fci,fc 2 ,fc,w) = r^-(fci,fc 2 ,fc,w) 

Nki.a ^ 3 k,k 2 Fa(k,(jj) 

V 5 (k 1 ,k 2 ,k,u]) = T^ a (ki,k 2 ,k,u>) 

~< ttfci.cr > &,fc 2 F a (k,Uj) 

T 6 (ki,k 2 ,k,u) = T^§ w (ki,k 2 ,k,uj) 

Tl kl ^a d k ,k 2 Ecr(k,id') 

(A5) 

which, by replacement in Eq. 1131 leads to: 




(fc,w) 


1 


2 vr [w - £ d (/c) - ]F fei (jj{N kl ,cr) + ^(nfci))] 

(A6) 


Finally, considering that: 


jy ^ 'A n kl,cr) ~ { n i,cr) — ( n c,cr ) Vi (A7) 

fel 

= w = <n d > w ( A8 ) 

/Cl 


where: n kl ,a =< >• 

Similarly, in this first level of approximation: 

r 2 (fci, fc 2 , k, lo) = T cdd a {k\ 1 k 2 , fc, w) 

> $k.k 2 a(k,Lu') 

T 3 (k 1 ,k 2 ,k,u) = T c a d 4-(ki,k 2 ,k,uj) 

~< N klilf > 5 k ,k 2 Ga{k,uj) (A2) 

where: -/V fcliCr =< <4 1)0 .<4 ll0 - >■ 

Therefore, replacing this approximation for Tj ( i = 1,3) 
in Eq. [T2]one has: 


one can rewrite the Green’s functions in this first order 
(Hartree-Fock) order of approximation as: 

G " F(fc,W) “ 2 t r[ U -E c (k)-(U(nc,a)-V(n d ))] 

(A9) 

and, 

^ ,UJ ^ 27t [oj - E d {k) - (U(nd,a) - V(n c ))] 

(A10) 


2. Second-order solution of EOM: RPA decoupling 


[uj - E c (k)]G a (k,uj) = ^■4 ll fc 2 (n fcl ,^)G CT (fc,w) 

ki,k2 

+ -j^dk 2 ,k ((N kl ,a) + (N kl ,a)) G a (k,U}) 

(A3) 

from which one directly obtains: 


Gf F (fc, W ) ~ ^ [w _ Ec{k) _ ^ ^ {nki _ } + F <JVfci) )] 

(A4) 

where: (JV fcl ) = (lV felt ) + (Nfcif)- 

On the other hand, using this approximation for ( i 


In the following, we will describe the second-order 
decoupling and solution of the equations of motion for 
G(fc, w) and F(k, w), which we used to calculate all spec¬ 
tral density and total density of states results presented 
in Section cm 

As mentioned above, the exact equation of motion 
for G(k,uj), Eq. [12] involves sums over two k-indices 
of the three coupled Green’s functions: Ti (i=l,3) de¬ 
fined in Eqs[l4] Now, instead of approximating these 
coupled Green’s functions directly in Eq. [T2] as in 
the first order (HF) solution, we determine the three 
equations of motion for the Tj (i=l,3) Green’s func¬ 
tions, respectively. Concretely, for r 1 (fci,fc 2 ,fc,w) = 

< Ck 2 ,<rc{ 1 ,cr c ki-k 2 +k,ir,c\ a > (w) we obtain the fol¬ 
lowing equation of motion: 
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u]T 1 (k 1 ,k 2 ,k,u;) = + 

< [Ti(ki,k 2 ,k),H\ ;4 ct > (w) (All) 


where pi, W] = pi, Ho] + Pi, V int \ introduces coupling 
of the EOM to new higher order Green’s functions. In 
particular, the term pi,"Ho] couples new Green’s func¬ 
tions << A;c» where A involves three operator 
products: one creation and two annihilation ones, e.g.: 

=< ^- C !' — k 2 -\-k,W c k ' 2 ,a'i 4,er ^ ( w )- 

On the other hand, the term pi, I/)nt] in Eq. IA11I cou¬ 
ples four new higher order Green’s functions, each in¬ 
volving five operator products in A: two creation and 
three annihilation ones, such as c a cf-Cacfca, c a clc a cf-Ca, 
Co-c^Cadl-da and Ca-cfrCad^-d-^. Similar results are ob¬ 
tained when calculating the equations of motion for 
T 2 (fci, fc 2 , fc, w) and r 3 (fci,fc 2 ,fc,w). 

Analogously, the exact equation of motion for F(k, to), 
Eq. 1131 involves sums over two k-indices of three coupled 
Green’s functions: Ti (i= 4,6) defined in Eqs[15l whose 


own equations of motion indirectly couple F(k, to) to new 
higher-order Green’s functions in the problem. 

To close the coupled set of equations of motion at this 
second-order level, we used the following approximation: 
all new higher-order Green’s functions introduced in each 
subset of equations of motion for T^ (i = 1,3), were 
approximated in mean-field in terms of (* = 1,3) and 
G(fc,w) (introducing appropriate average values), yiel¬ 
ding a 4 x 4 closed set of equations of motion. We pro¬ 
ceeded likewise for the subset of equations for (i = 4,6) 


related to F(k,co). 

To exemplify the decoupling scheme adopted to close 
the set of EOM at this level, below we show the approxi¬ 
mation used for the higher-order Green’s functions in¬ 
volving three operator products: k 2 , k, oj) = 

< Ck' 2 ,ac{' itW c k ^ k2+kjW ] c' k (T > (w), according to the 
definition in Eo llfil Approximating this Green’s function 
in mean field yields: 


y\cdd 

1 CT,<7,C7 


{k[,k 2 ,k,uj) 


(4i,a C fci-fc2+fc,a) < C k' 2 ,cr'i C fc,cr > M 
d k / 2 k P/q,7f) C k ' 2 ,cn C k(7 (oj) 

,7t)G a ( k , Ca) 

I 


Equation of motion for Ti(fci, k 2l k, w): 


We now illustrate the decoupling of the 
higher-order Green’s functions involving 
five operator products, by two examples: 

c k' 2 ,a c k i it -ff C k‘i-k' 2 + k ' 1 ,aC^ 1 -ffC kl — k2 -^ k ^\C ka 

and <!C c k ' 2 ,<T c k ^^ c k 1 —k 2 +k,'ad k ^ (J d k2 — k ' 2 ^ k ' i ^ a ; c ka i§>. 

Approximating these Green’s functions in mean field 
yields: 


^ C k' 2 ,a c k ' i ^ c k2-k , 2 +k' 1 ,wCl l W Ck 1 -k 2 +k,a', 4 ,ct ( UJ ) 

— ( c k' 1 ,w c k2-k 2 +k' 1 ,a) ^ c k' 2 ,(r c kl ^ c ki—k 2 +k,a i c ka ^ ( w ) 
”l _ ( C ii,CT Cfe i-fc 2 + fc, 5 : ) C k ' 2 ! crC ki ~C kl — k 2+k,W i c fc,er --® > ( W ) 

- s kl>k > a {rik^T^^kuk^k, w) 

+5 kM {n k ^)T^{k[XM 

— 5k 1 ,k' 2 {nk’ 1 ,o : )ri{k 1 ,k l 2 ,k,u]) + S k ,k 2 (n kl ,w)^i{k[,k 2 , k,u>) 


and, 


^ c k 2 ,a c k[,a C kl—k 2 +k,ad k ^ a d k2 — k ' 2 + k ^ a ; C fe a ( w ) 

— { C \ 1 ,a Ck l-k 2 + k,a) C k ' 2 ,<r< 4 ' )(T ^fc 2 -fc'+fc',<r i c k,a ^ ( W ) 

^'^k' 1 ,a < ^k2-k 2 +k , 1 ,a) ^ c k 2 ,ir c kl t a C ki — k 2 +k,W i c fc,<7 (w) 

— d k2 , k (n kl ,tj)T 2 ( k 3 , k 2 , kj ca) -f- S k2 k ' 2 (Aq.* j( 7 )r 3 (& 3 , k 2 , k : cc) 


All other higher-order Green’s functions were decou¬ 
pled following the same procedure. We proceeded like¬ 
wise for the subset of equations for Tj (i = 4,6) related 
to F(k, to). 

Finally, the equations of motion for the six Green’s 
functions Tj (i = 1,6) in RPA approximation read: 


wTi(ki,k 2 ,k,w) = —^—(n k i,w) + [E c (ki -k 2 + k) - E c (k\) + E c (k 2 ) + U{n c + l)]V 1 (ki,k 2 ,k,uj) 

Z7T 

{21^ {7lki,a)^ld H - (2 + (jlk-L Tl c + U {j^k\ — k2+k,w)) }^cr (&? 

y^r 2 (fc 1 ,fc 2 ,fc,o;) + y^r 3 (fc 1 ,fc 2 ,fc,a;) 

m k 1 k 1 


V 

+ ~jy ({n-ki,w) — ( klk 1 -k 2 +k,a )) 


(A12) 


Equation of motion for r 2 (£q, k 2 , k,u >): 
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uT 2 (k u k 2 ,k,u) = + [E d ik\ - k 2 +k)~ E d (ki) + E c (k 2 ) + Un c + 2n d (U + V) - V]r 2 (Aq, k 2 ,k,ui) 

z 7 T 

+{([/ - V)(N kua ) + 2 V(N kua )n d - V((N kli<r ) - (N kl _ k2+k , a )) + V}G a {k,u) 


+ jy {{E kl ,rr) {Nk 1 -k 2 +k,a)) 


^2r 1 (k 1 ,k 2 ,k,w) 


+ jy i(Nki,cr} {N kl — fe 2 +fc,cr)) 


^ r 3 (fci,fc 2 ,fc,w) 


(A13) 


Equation of motion for T 3 (kpk 2 ,k,uj): 


u>T 3 (ki,k 2 ,k, uj) = ^p^{N klW ) + [E d (ki - k 2 + k) - E d {k\) + E c (k 2 ) + Un c + 2 Vn d + V]V 3 (ki,k 2 ,k,uj) 

Z 7 T 


+ {U{N kl ,cr)n c + G (1 {ri kl ^ a ') (j'l kl - k 2 -\- k .( 7 ))'IT d T C < j ) 


V 


+ ((Afc llC r) (Afci-fe+fejO-)) 


r i(fci, A: 2 , A;, oj) 


+ {(Nki,a) {Nk 1 -k 2 +k,a}) 


'^2r 2 (k 1 ,k 2 ,k,u) 


(A14) 


Equation of motion for r 4 (/q, fc 2 , fc,w): 

wr 4 (/ci,/c 2 ,fc,w) = %^-(AT fc *) + [£ d (fci - k 2 + k) - E d (k i) + E d (fc 2 ) + J 7 (n d + 1 )] r 4 (Aq, k 2 , k, uj) 

Z 7 T 

+ {2V(N kl , w )n c + (2 + ( N kl , w ))Un d + U((N kli5f > - (JV fcl _ fca+fci? ))}F (r (fc,u;) 


+ ^ {(N kl ,<r) (N kl — k2 + kj(T )) 


TM^M + YMk^M 


(A15) 


Equation of motion for ^(Aq, k 2 , fc, uj): 


LdT 5 (ki,k 2 ,k,cv) = ^ 7 p^{n kltl7 ) + [£? c (fci - k 2 + k) - E c (ki) + E d (k 2 ) + Un d + 2n c (U + V) - V] T 5 (k 1 ,k 2 , k,u) 

Z 7 T 


+{(f/ - V)(n kl , a ) + 2 V(n kl ,a)n c - V{{n ku<J ) - ( n kl - k2+ktCr )) + V}F CT (A;,w) 


V 


+ ^ ({ n ki,cr) { n ki~k 2 +k,cr)) 


r 4 (fc 1; fc 2 , k, uj) 


+ ^ ({ n ki,cr) ( nki-k 2 +k,cr }) 


'^2r 6 (k 1 ,k 2 ,k,uj) 


(A16) 


Equation of motion for Ve(ki,k 2 ,k,uj): 


r 

wrg(fc 4 , k 2 ,k,u j) = 2 ’ (rtfc! j) + [F c (Aq — k 2 + k) — E c (k{) + E d {k 2 ) + Un d + 2Vn c + V] r 6 (fci, k 2l k, uj) 

Z7T 

+{U(n kl ^)n d + V{1 - {N kl ^) - (. N kl _ k2+ktW ))n c + V((N kl ^)}F a (k,uj) 


V 

+ "^ {{ n ki,a) ~ { n ki—k 2 +k,a}) 


^2r 4 (k 1 ,k 2 ,k,u}) 

k‘ 


+ (( n ki,(T ) ( n ki—k 2 +k,tr}) 


i,k 2 ,k,u) 

k’ 


(A17) 


Notice that this RPA-like approximation involves a second-order closed set of equations for G a (k,uj) 1 and double 
k-index summations of Id, T 2 and Tg; and an analogous one for F a (k : uj), T 4 , Ts and IV Below we explain how 
we could solve these two sets of coupled equations of motion, and detail the Green’s functions G(k,uj) and F(k,w) 
obtained, required to calculate the spectral density and total density of states results presented in Section [ 


a. Solution of the equations of motion determining G a (k, to). 

The keypoint to close the EOM system at this second-order level, enabling to determine G a (k,uj ) from Eq. [T21 
using the RPA-coupled set of Eqs. IA12IA13l an d lA141 is that only BZ summations over Aq and k 2 of the P,;(Aq, k 2 , k) 
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(i=l,3) functions appear in Eq. [12] So that, appropriately (double) summing over the BZ the set of coupled set of 
equations, one can effectively extract G a (k,uj) at a fixed k as we explain in the following. 

It is useful to introduce the following notation: 


<1 >j(k2,k,uj) = y^T i (k 1 ,k2,k,ui) , * = 1,3 

ki 


which allows to rewrite Ecis lA12l IA13I and IA14I as follows: 

r 1 (k 1 ,k 2 ,k,u) = A\ + BiG a ( k , uj) + G\ 
r 2 (fci, fc 2 , fc, uj) = A 2 + B 2 G a(k, uj) + C 2 

r 3 (fci, k 2 , k,u>) = A 3 + B 3 G a (k, to) + C 3 


[$2 (fc 2 , k, uj) + $ 3 {k 2 , k, w)] 


where the coefficients A,, Uj and Cj are: 


<h 3 (fc 2 ,fc,w) 

$ 2 (fc 2 ,fc,w) 


A 1 


A 2 

A 3 


(**fci ,<r) 

27t [w — E c (k 2 ) — t/(n c + 1)] 

_(iW)_ 

2vr [w - B c (k 2 ) - Un c - 2n d {U + V) + V] 

_(jV fcl , g )_ 

27 t [uj — E c (k 2 ) — Un c — 2Vrid — V”] 


(A18) 


(A19) 


73 { 2 E(rifej } (j)Tid (2 A ,a))Un c A £/((n/ Cli( y) 

— - 

W — Wi 

„ {(E/ - V)(N kl , a ) + 2 V(N kl ,cr)n d - E«JV fcli(7 > - (Af fcl _ fc2+fe , CT )) + E} 

i> 2 = - 

W — W 2 

{{/ (N klt „)n c + E(1 - (n kl ,a) ~ {nk 1 -k 2 +k,a))n d + V((n kl ,o’)} 

B 3 — - 

UJ — UJ3 


and, 


Ci 

c 2 

C 3 


N [( n ki,a) { n k\ — k2 + k,a)} 

UJ — UJ\ 

l(N kl ,cr) — (N^-kz+k^)] 

UJ — UJ 2 

[(A^fci.g) ~ (N kl -k2+k.a)} 

UJ — UJ3 


uj 1 = [ E c (k\ — fc 2 + k) — E c (k\) + i? c (fc 2 ) + U(n c + 1)] 

uj 2 = [Ed(k\ — fc 2 + k) — E d {ki) + E c (k 2 ) + En c + 2n d {U + E) — E] 

uj 3 = [£'d(^i — fc 2 + k) — Ed{k\) + E c (k 2 ) + t/n c + 2E n d + E] 

Notice that if Eas lA19l are BZ summed over k\, each of the <hj(fc 2 ,fc) can be derived in terms of G a (k,uj) and 

coefficients involving the temperature dependant-fillings, the effective band structure and the model’s correlation 
parameters. 

Meanwhile, if one introduces the notation: 

Qi(fc) = Y. Tj(ki,k2,k,Qj) = YMk^k, uj) , * = 1,3 (A20) 

fcl,/C2 Al2 

Eq. QJ can be rewritten as: 


1 

27T 


^0i (k,u) + ^e 2 (fc,w) + ^e 3 (k,uj) 


[uj - E c (k)] G a (k,uj) 


(A21) 
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Performing the k 2 summations in Eq. IA201 all @i(fc) can be expressed in terms of G a (k , w) and coefficients involving 
the temperature dependant-fillings, the effective and the model’s correlation parameters, so that one can finally solve 
for G a (k,uj). We find the following expression for G (k, w), related to the c-electrons in our model, in this RPA- level 
approximation: 


G^ A (k,u)^ 


h + (xr) + 7r it + [^- 3 ]fci + [G 3 ]fci (jfc + jt) } 


w -E c {k) - Y,k 2 {77 (fr) + F %. + t-^k + [C 3 ]fci (f^- + If) } 


where, denoting: [A]*^ = A, one has: 


(A22) 


= ( 1 - ^[CfekPak) (1 - [CilfcJCak) - [C'lJjfcx (^k + l) (^[Cak + [Cak 

Yi(fc 2 ,fc) = ( 1 — jy[Cak[C , 3]fe 1 ^ ([Ai]fcj + [A 3 ]fe 1 [CiJfcJ + ^[A 2 k + ■^[C , 2]fc 1 [Cik ^jyPak + 1 

Zi(k 2 ,k) = ( 1 - -^[C 2 ]fc 1 [C' 3 ]fc 1 j {[Bi] kl + [B 3 ] fcl [Ci] fcl ) + [-B 2 k + ■jy[C' 2 ]/ 5l [B 3 ]fc 1 


*a(fca,fc) = (1 - [CrkPak) ( 1 ' ^[CakPsk) - (jj[C 3 ] kl + V ) [C2) kl [Ci\ kl ({C 3 ] kl + 1) 

Y 2 {k 2 ,k) = (1 - [Cik[C 3 k) ( [A*]* + jylCakl/kkI + ^ —[C 3 ] fcl +V ) [C 2 ] fel ([Ai] fcl + [A 3 ] fcl [Cik) 
Z 2 (k 2 ,k) = (1 - [CikPak) ( [B 2 ] kl + ^[C 2 ] kl [B 3 } k A + (^[C 3 ] kl +V] [C 2 ] kl + [B 3 ] kl [C 1 ] kl ) 


b. Solution of the equations of motion determining F tT (k, uo). 
We proceeded likewise for the subset of equations for Ti (i = 4... 6 ) coupled to F a (k,uj): 

r 4 (fci, fc 2 , k,uS) = A* + B*F a (kiui) + C* [d> 5 (fc 2 , k, w) + <h 6 (fc 2 , fc,w)] 

r 5 (fci, fc 2 , fc, w) = A 2 + w) + C 2 * 

r 6 (fcr, fc 2 , fc, w) = A 3 + w) + C 3 * 


jA $ 4 (k 2 ,k,u) + f^A d> 6 (fc 2 , k,w) 


jA $ 4(fc 2 ,fc,o;) + (^A $5(k 2 ,k,uj) 


(A23) 


where, 


AJ = 
A* = 


(N kl ,a) 


27r [w - E d (k 2 ) -U(n d + 1)] 

{tiki ,cr) 

2 = 2tt [w - E d {k 2 ) - Un d - 2n c {U + V) + V] 

A* _ _ ( n ki ,0) _ 

3 “ 2tt [w - E d {k 2 ) - Un d - 2 Vn c - V] 


B{ = 


B* = 


b; = 


{2V(N kl ,y)n c + (2 + (N kl ,g))Un d + U{(N kl , g ) - (JV fel _ fc2+fc , g ))} 

U! — U>1 

{{U - V){N kl , a ) + 2 V(n kl , a )n c - V((n ku<7 ) - (n kl - k2+k , a )) + V} 

ui — w 2 

{U(n kl ,w)n d + V(l- (N kl ,w) - (N kl _ k2+k ^))n c + V((N kl ,w)} 

W — Wo 
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C* 

C* 2 

c; 


V [(-^fci.g) - (Nk 1 -k 2 +k,a)] 
N uj — oji 

[( n kl ,a) ~ (^fci-fc 2 +fc,cr)] 

UJ — UJ 2 

[(n kl ,w) ~ ( r lk 1 -k 2 +k,a )] 

UJ — UJ 3 


w* = [E d (ki - k 2 + k) - E d {k\) + E d {k 2 ) + U(n d + 1)] 

uj 2 = [F c {k\ — k 2 + k) — E c (k \) + E d (k 2 ) + Un d + 2n c (U + V) — V] 

Wg = [S c (fci - k 2 + k) — E c (ki ) + E d (k 2 ) + Un d + 2Vn c + V} 

Carrying out the summations over k\ and k 2 BZ points in the set of equations in Eq. IA23l and Eafl3l we analogously 
find an expression for the Green’s function F a {k 1 to) corresponding to the d-electrons, in this RPA-level approximation: 


F RPA (k,uj) = 


f U 

(y;> 

1 ^ r 

\ N 


+ N 


XL 

x* 


+ 


i A 3.]ki + [Oik + Jc![) } 


u 

(z r> 


Ik* 

N 


+ N 

X 2 


+ 


[-®3 ]ki + [<^3]fci (x* + X*)] } 


(A24) 


where: 


Xt(k 2 ,k) = ( 1 - ^[c* 2 ] kl [c* 3 ] kl ) (1 - [c?k[<? 3 *k) - [c?k 

u, 


N l 

U, 


Fi (k 2 ,k) — ( 1 AT [C 2 ] kl [G 3 \ kl ) ([kk + [^ 3 ] fci [Ci k) + ( kk + AT [C 2 \ kl [kk ) ( ^[Cgk + 1 


[Cg]fci +1 


U, 


U V 

N^ + N 


N l 

U, 


[Calfci 
U , 


Zl(k 2 ,k) = ( 1 - jj[CZ} kl [C* 3 ] kl ) (kk + kkkk) + kk + - T [C* 2 \ kl [B;] kl 


N l 


x; (k 2 ,k) 


Y 2 *(k 2 ,k) 


Z* 2 (k 2 ,k) 


(1-kkkk) 
(1-kkkk) 
(1 -kkkk) 


(l - ^klkksk) - (^kk + V^j [CS\ kl [CtW «C 3 *k + 1) 

(kk + |k 2 *kkk) + (Jj[c* 3 ] kl + V^j [c 2 *k (kk + kkkk) 

(kk + ^kkkk) + (^kk + V s ) [C* 2 ] kl (kk + kkkk) 


Appendix B: Analysis of the k-dependent self-energy 

To better assess the nature of the second order (RPA- 
like) approximation which we used to decouple the equa¬ 
tions of motion and determine the relevant Green’s func¬ 
tions, yielding the spectral density and total density of 
states of our model, in the following we present an ana¬ 
lysis of the self-energies for c and d electrons obtained 
in our approach, focusing on the k-dependence and the 
temperature dependence. 

Using Dyson’s equation, one can determine the second 
order self-energy related to each electron band as follows: 

E c (jfe, u) = G°(k, uj)- 1 - G RPA (k , w )" 1 (Bl) 


(B2) 


where, G°(fc,w) and F°(k,oj) represent the retarded non¬ 
interacting Green’s functions corresponding to the two 
bare effective bands: for the uncorrelated c- and d-bands, 
respectively, G°(k,uj )~ 1 = u — E c (k) and F°(k,uj)~ 1 = 
u> — E d (k). Meanwhile: G RPA (k, uj) and F RPA (k, uj) were 
defined in Eqs. A22 and A24 of Appendix A. 

First, notice that the self-energy in Hartree-Fock 
approximation for the c- and d-orbital Green’s functions 
(Eqs. A4 and A 6 ) is independent of crystal momentum 
k and energy. 

Therefore, in this appendix we present our results ob¬ 
tained evaluating in second order level of approximation 
(as detailed in Appendix 0 the real and imaginary parts 
of the self-energies defined above, corresponding to the 
retarded electron Green’s functions of eqs. A22 and A24. 
Our results confirm that in our RPA-like approach, the 
self-energies exhibit clear renormalization effects due to 


S d (k,u) = F°(k,uj)~ 1 - F RPA {k,uj)- 1 
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the correlations included: with non-trivial crystal mo¬ 
mentum, temperature, doping, and energy dependence. 



ffl(eV) 


FIG. 8. Imaginary part of the c- and d-band self-energies as a 
function of energy (w.r. to the Fermi level), at the following 
symmetry points of the square lattice Brillouin zone: r ( fc = 
(0,0)), M{ k = ( 7 r, 7 r)) and X ( k = ( 7 r, 0 )). Parameters: 
U = V = 3.5, n c = 1.45 and rid = 0.56, T = 20A'. Bare 
effective band tight-binding parameters as in Ref. [3^ 


In the following, our results will be presented in a se¬ 
ries of figures, in particular exhibiting the diversity of 
the temperature dependence which is obtained at diffe¬ 
rent k-points in the Brillouin zone. As general consis¬ 
tency checks of our results, notice: the correct (neg¬ 
ative) sign exhibited by the imaginary parts of the 
electron self-energies, while the real parts of the self¬ 
energies, which represent the energy shifts renormalizing 
the bare electron energies in our approximation, are cor¬ 
rectly Kramers-Kronig related to the respective imagi¬ 
nary parts. 

1. Momentum dependence of E c (fc,aj) and S d(k,u)) 

First, we exemplify the momentum dependence of the 
imaginary part of the two self-energies E c (k,w) and 
Hd(k,ui) in our approximation. In Fig. [51 we show the 
imaginary parts of the two self-energies at the three rele¬ 
vant high symmetry BZ points of the square lattice, usu¬ 
ally probed by ARPES. The upper two panels show that 
at the Brillouin zone centre (r) and at M, the c-band 
is mainly renormalized around the Fermi level, while the 
d-band is mostly renormalized far away from the Fermi 
level, at energies close to the bare band edges (w > 5eV). 


The lower panel depicts the imaginary parts of the two 
self-energies at X (k = (-7T, 0)): which noticeably differs 
from the previous cases. At X the renormalization is re¬ 
levant mostly near the band edges of the bare electronic 
structure, while near the Fermi level the renormalization 
at X is not significant. 

In the rest of this appendix, we will focus on the renor¬ 
malization of the energies obtained near the Fermi level, 
and therefore center our discussion on E c (fc,w). 



FIG. 9. Temperature dependence of the imaginary part of 
E c (fc,u;) at T. T as labelled in the figure, other parameters 
as in Fig. [ 8 ] 


2. Temperature dependence of the self-energy 

a. Imaginary part of the self-energy 

Now, we will focus on the evolution with temperature 
of the imaginary part of E c (fc,w), in particular fixing k 
at the relevant high-symmetry points of the BZ studied 
by ARPES experiments, and also at other BZ points: 
in which we found that the temperature dependence of 
the renormalization was larger and non-trivial (not the 
expected behaviour resulting from /j?c(u)). 

Fig. [9] shows that at T the renormalization of the self¬ 
energy is almost independent of temperature. We found 
that also at M and X the effect of temperature on the 
renormalization is almost irrelevant. This negligible tem¬ 
perature dependence of the self-energy at the high sym¬ 
metry BZ points explored by ARPES, T, M and X , ex¬ 
plains the spectral density function results presented in 
Figs. 10 and 11 of Section ??. 

But from the total density of states results in Fig. 
9 of Section ??, we know that non-trivial relevant 
temperature-dependent effects indeed appear when the 
contributions of the whole BZ are taken into account. 
Therefore, we used a grid to explore the temperature 
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FIG. 10. Temperature dependence of the imaginary part of 
E c (fc,w) at the BZ points respectively indicated in the plots. 
T as labelled in the figure, other parameters as in Fig. [ 8 j 


FIG. 11. Temperature dependence of the real part of E c (k, ui) 
at the BZ points respectively indicated in the plots. T as 
labelled in the figure, other parameters as in Fig. [ 8 ] 


dependence of the renormalization at different points of 
the BZ (528 points), and interestingly we could identify a 
number of specific BZ points, not yet probed by ARPES, 
where sizeable non-trivial temperature dependent renor¬ 
malization effects are obtained. 

Concretely, in Figure llOf al we exhibit the tempera¬ 
ture dependence of the imaginary part of E c (fc,w) at 


k = (0.557T, 0.557 t): a BZ point where we find that, 
increasing temperature, interaction-related renormaliza¬ 
tion effects lead to a redistribution of spectral weight 
near the Fermi level, with temperature dependent peaks 
evolving in the energy range [—0.4, —0.2] eV. 

In Figure flOl' bl. we exhibit the temperature de¬ 
pendence of the imaginary part of £ c (fc, w) at k = 
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(0.67T, 0.67r), where the spectral weight redistribution is 
even larger, involving e.g. a number of states present 
above the Fermi level at T= 40K , which are pushed be¬ 
low the Fermi level at higher temperatures, leading to 
a number of peaks with temperature-dependent location 
and height in the energy range [—0.7,0] eV. 

Finally, in Figure me) we exhibit the temperature 
dependence of the imaginary part of H c (k,ui) at k = 
(0.567T, 0.3l7r). Notice that, increasing temperature, the 
interaction effects lead to temperature dependent peaks 
now in the energy range [—0.6,—0.1] eV. It would cer¬ 
tainly be interesting to investigate these predictions by 
ARPES. 


b. Real part of the self-energy 

To complete the presentation of the renormalization 
effects described in our approach, in Figure [IT] we 
show the evolution with temperature of the real part of 
E c (£:,w) at the BZ centre, where it is negligible, and at 
two of the k-points we identified as having relevant non¬ 
trivial temperature dependence, discussed above. The 
results in Fig. [FT] are the Kramers-Kronig related coun¬ 
terparts of the imaginary part of the self-energies shown 
in Figs. [9] flOlal and HUT bl respectively. We checked 
that the renormalizations at M and A', as expected, are 
almost independent of temperature. 


Appendix C: Description of doping effects on the 
electronic structure 

Here, we complement our discussion of the effects of 
doping on the electronic structure of Section 3.3, mainly 
presenting spectral density function results. In particu¬ 
lar, we analize the main changes in A(k 7 to) with doping, 
fixing k at the relevant high-symmetry points of the BZ 
which were experimentally studied. We compare with 
ARPES [76] and angle integrated valence band photoe¬ 
mission experiments [77] on undoped, Co-doped (electron 
doped), and K-doped (hole doped) Ba-122 single crystals, 
in which asymmetric electron-hole doping effects on the 
electronic structure were reported, [76,77] confirming pre¬ 
vious indications of Hall and transport experiments. [79] 

In Fig. H2J we show A(k,uj) for three cases, corre¬ 
sponding to the fillings of the Ba-122 compounds stud¬ 
ied in Refs.[76,77]: the undoped (n = 2), hole-doped 
(n = 1.85), and electron-doped (n = 2.25) systems, at 
the r and M Brillouin zone points. As observed at 
r with ARPES,[77] we obtain two main peaks (apart 
from a few other peaks all of lower intensity): one of 
them centered at binding energy 0.55 eV for the un¬ 
doped compound, and the other one near the Fermi level, 
which shift their positions according to doping follow¬ 
ing the trend observed in experiments. Concretely, the 
peak at u> ~ —0.55 eV for the undoped system, appears 



-0,8 -0,6 -0,4 -0,2 0 0,2-0,8 -0,6 -0,4 -0,2 0 0,2 

o) (eV) 


FIG. 12. Comparison of the effect of electron vs. hole doping 
on A(k,uj), at T and M points. T = 40 K, U = V = 3.50eV, 
v = 9. Other parameters as in Fig.l. 



co ( eV ) 


FIG. 13. Effect of a) hole doping and b) electron doping on 
A(k,cj), at M point. T = 40 K, U = V = 3.50 eV. v = 9. 
Other parameters as in Fig.l. 


shifted away from the Fermi level: at u ~ —0.60 eV 
for the electron-doped system, while for the hole-doped 
system it appears shifted in the opposite direction: at 
u) ~ —0.50 eV. In fact, at T the whole spectrum ap¬ 
pears shifted analogously as a function of doping. Thus, 
electron-hole asymmetric effects such as experimentally 
observed[76,77] are present in the renormalized electronic 
structure of our model. Regarding the spectral densities 
at M shown in Fig. [121 the situation is seen to differ: 
though the same trend of spectral weight shifts upon 
doping sign is seen for the part of the spectrum corre¬ 
sponding to energies closer to the Fermi level (oj > —0.4 
eV), an important spectral weight redistribution is also 
present upon doping. 
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FIG. 14. Total DOS A(w), for different values of hole doping 
(indicated in the figure) at T = 20 K, 17 = V = 3.50eI7. Other 
parameters as in Fig.l. Inset: dependence of the chemical 
potential on total band filling. 


In Fig. [13] we focus on the spectral density at M in 
the range of energies closer to the Fermi level studied 
by ARPES in Ref. [76], and exhibit our results for more 
values of charge doping: (a) for hole doping, and in (b) 
for electron-doping. The nature of the spectral weight 
redistribution that we find takes place upon doping is 
in qualitative agreement with the reported data. [76] In¬ 
creasing hole-doping, spectral weight from lower energies 
is transferred towards the Fermi level. Meanwhile, in¬ 
creasing electron doping we obtain a smooth increase of 
the quasiparticle weight at the Fermi level, while other 
states are pushed away to lower energies. 

In Figure 1X71 we complement the results presented in 
Fig. 3, showing the total DOS for a larger number of 
hole doping values One clearly appreciates that increas¬ 
ing hole-doping, a rigid-band-like shift of all peaks to¬ 
wards the Fermi level is obtained, thus increasing spec¬ 
tral weight there in agreement with ARPES results on 
doped BaFe2As2,[76] where a smooth chemical poten¬ 
tial shift upon doping, related to the effective masses of 
the low-energy valence states, was observed. The filling- 
dependence of the chemical potential in our model is 
shown as an inset in Fig. [14] 

Finally, in Figure C.8 we exhibit effects of doping on 
the renormalized Fermi surface (FS) topology of the cor¬ 
related two-orbital model. Notice that several Lifshitz 
transitions are obtained. Starting from the FS typical of 
parent compounds (n = 2), doping with electrons mainly 


reduces the hole point around the BZ centre T, until it 
disappears at critical electron filling n e ~ 2.65. At the 
same time, one can see that electron doping increases 
the four electron pockets located around the BZ points: 
X,Y, -X,-Y, until a new FS topology emerges at n e , with 
a larger electron FS around T. This is qualitatively 
consistent with experimental reports [22,23]. Similarly, 
notice in Figure C.8 that hole doping reduces the four 
above-mentioned electron pockets until at critical hole 
filling rih ~ 1.81 they disappear, and a new FS topology 
emerges: only consisting of the hole pockets. 

1. Temperature dependence of the electronic 
structure: effect of doping. 

After having discussed in section 3.4 the effects of tem¬ 
perature on the electronic structure for the parent com¬ 
pounds ( n = 2 ), here we compare them with the results 
predicted for doped systems. 

In Figure C.9 we focus on the total density of states: 
the results shown evidentiate that the non-trivial tempe¬ 
rature dependent renormalization effects by correlations 
found for the parent compounds, shown in Fig. C.9.(b) 
and previously discussed in connection with Figure 4, are 
greatly reduced when doping is introduced. In particular, 
Fig.C.9.(a) shows that for hole doping, increasing tem¬ 
perature, the dominant peak near the Fermi level is not 
displaced and the main effect of temperature is the usual 
expected thermal broadening (and consequent decrease 
of the height of the peak). While for electron-doping, 
Fig.C.9.(c) shows that a certain non-trivial redistribu¬ 
tion of spectral weight still appears with temperature, 
though much less important than for the parent com¬ 
pounds: notice that the dominant peak near the Fermi 
level, is shifted to lower energies by an increase of tem¬ 
perature, reducing its height. 

We also explored the effect of doping on the tempe¬ 
rature dependent spectral density function results. In 
Figure C.10, for an electron-doped system, we focus on 
the spectral density function: in particular for the three 
Brillouin zone points previously discussed in connection 
with Figure 6, exemplifying BZ points we could iden¬ 
tify for the parent compounds where larger non-trivial 
temperature dependent renormalization effects by corre¬ 
lations are predicted. Notice that doping strongly affects 
those k-dependent results: for the two BZ points shown 
in Figs. C.10.(a) and C.10.(c), the presence of doping 
has suppressed the temperature dependence. While Fig. 
C.10.(b) shows that for &2 = (0.67T, 0.67t) some tempe¬ 
rature dependence is still retained. For hole doping, we 
found that no relevant temperature dependent effects are 
retained at these three BZ points. 
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FIG. 15. Correlated two-orbital model: doping evolution of the (large) Fermi surface (one Fe atom / cell), at T = 40 K. Electron- 
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FIG. 16. Total DOS A(w) shown at three temperatures: T = 
10 K, 60 K and 100 K, for fillings: (a) n = 1.75 corresponding 
to hole-doping; (b) n = 2 ( parent compounds); (c) n = 2.25 
corresponding to electron-doping. U = V = 3.5eV, v = 8. 
Other parameters as in Fig.l. 


FIG. 17. Temperature dependence of A(k,ui) for n = 2.25: 
shown at three BZ points, respectively indicated in each plot. 
Parameters: U = V = 3.5 eV , v = 8. Other parameters as 
in Fig.l. 
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